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Abstract 

We study non-Gaussianity generated by adiabatic and isocurva- 
ture primordial perturbations. We first obtain, in a very general set- 
ting, the non-linear perturbations, up to third order, for an arbitrary 
number of cosmological fluids, going through one or several decay 
transitions. We then apply this formalism to the mixed curvaton 
and inflaton model, allowing for several decay channels. We compute 
the various contributions to the bispectrum and trispectrum resulting 
from adiabatic and isocurvature perturbations, which are correlated 
in general. By investigating some hybrid decay scenario, we show 
that significant non-Gaussianity of adiabatic and isocurvature types 
can be generated without conflicting with the present isocurvature 
constraints from the power spectrum. In particular, we find cases 
where non-Gaussianity of isocurvature origin can dominate its adia- 
batic counterpart, both in the bispectrum and in the trispectrum. 



1 Introduction 



Observations of cosmic density fluctuations, such as cosmic microwave back- 
ground (CMB) and large scale structure, are our only accessible windows 
onto phenomena that occured in the early Universe. It is therefore crucial 
to extract as much information as possible from these observations, as a way 
to test or constrain various early universe scenarios (see e.g. [1] for recent 
lecture notes). 

This is the reason why primordial non-Gaussianity has been the subject of 
an intense study recently, as it would provide invaluable information beyond 
the power spectrum. In particular, a detection of primordial non-Gaussianity 
would rule out the simplest models of inflation (with a single field), which 
generate only very small non-Gaussianity. 

Observations of the CMB anisotropies by the WMAP satellite [2] have 
set the present limit f$ft? = 32 ± 21 (68% CL) [or -10 < f§? < 74 (95% 
CL)] on the non-linearity parameter /nl^ which characterizes the amplitude 
of the bispectrum, for the local shape* 1 . Although Gaussian fluctuations are 
still consistent with current observations, this constraint suggests that pri- 
mordial fluctuations may deviate from Gaussianity since the central value 
of this range is somewhat away from zero. The Planck satellite should be 
able to confirm or infirm such level of non-Gaussianity. If Planck detects sig- 
nificant local non-Gaussianity, then early Universe scenarios with additional 
fields, such as a curvaton [3] or a modulaton [4-9], would become natural 
alternatives to the simplest inflation scenarios. 

As soon as one considers models with multiple scalar fields (several infla- 
tons or spectator fields beside the inflaton), one must take into account the 
possibility that some (baryon or dark matter) isocurvature perturbations may 
be generated, in addition to the usual adiabatic fluctuations. The amplitude 
of such isocurvature modes, which could be partially or fully correlated with 
the adiabatic one, is now severely constrained by observations of the CMB 
power spectrum. However, complementary information, or constraints, on 
isocurvature modes can be obtained from non-Gaussianity. 

Non-Gaussianity in models with isocurvature fluctuations has been stud- 
ied in several works [10-15]. In particular, a general treatment for the primor- 
dial bispectrum which allows various decay scenarios has been given in [15] 
recently. Although the works mentioned above have mainly studied the bis- 
pectrum (except [14]), the trispectrum will also become, in the near future, 
a major target since combined information from the bispectrum and the 

# 1 Three types of /nl have been discussed in the literatures. Other two types and their 
constraints are —214 < < 266 for the equilateral type and —410 < /^ r L thog < 6 for 

the orthogonal type (95% CL.) [2]. 
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trispectrum is a powerful way to discriminate between scenarios that gen- 
erate local non-Gaussianity, as discussed in detail in [16]. In light of this, 
a general treatment of non-Gaussianity in models with isocurvature fluctua- 
tions for the trispectrum is worth investigating and is the main topic of this 
paper. More specifically, in the present work, we extend the general formal- 
ism given in [15] to third order, then apply it to the mixed curvaton and 
inflaton scenario [17-21]. 

The organization of this paper is as follows. In the next section, we present 
the general formalism to describe the curvature and isocurvature perturba- 
tions, up to third order. Then, in Section 3, we apply this formalism to the 
mixed curvaton inflaton scenarios and compute the bispectra and trispectra 
generated by the adiabatic and isocurvature modes in these scenarios. Some 
quantitative discussions are also given. The final section is devoted to the 
conclusion. Throughout this paper, we set the reduced Planck energy scale 
Mpi to be unity. 



2 General formalism 

In this section, we review and extend the general formalism introduced in 
[15] to describe the adiabatic and isocurvature perturbations generated in 
scenarios with a cosmological transition due to the decay of some species of 
particles, which we call a. Explicit results up to second order were given 
in [15], whereas we extend here the analysis up third order, in order to be 
able to compute the relevant trispectra. 

The purpose of the present analysis is to determine, at the non-linear 
level, the perturbations of an arbitrary number of fluids, after a cosmological 
transition where one of the fluid, labelled a, decays. For any fluid A, one 
can define, in a covariant way [22,23] (see also the recent review [24]), a fully 
non-linear curvature perturbation covector as 

where pa and Pa are energy density and pressure for the fluid A and Af = 
/ drV is the number of e- folds along the fluid worldlines (with proper 
time r) with being the four-velocity of the fluid. 

If wa = Pa/ Pa is constant, which will be assumed in the following as 
we will later consider only relativistic species (w = 1/3) and nonrelativistic 
species (w = 0), the covector (1) can be written as the total gradient of 

( A = 5AT+ ] ( 2 ) 
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where 5Af = Af — Af and the barred quantities are defined in a reference 
spacetime, which is strictly homogeneous and isotropic. The fully non-linear 
perturbation (a, which coincides with the definition of [25], is conserved on 
large scales for a non-interacting fluid. 

We also introduce the non-linear isocurvature perturbation of any fluid 
A, with respect to the radiation fluid, as 

S A = 3(( A -C r ), (3) 

where ( r is the curvature perturbation of the radiation fluid. 

Let us now focus on the transition due to the decay. In the sudden 
decay approximation, which we adopt here, the decay takes place on the 
hypersurface characterized by the condition 

H d = r CT , (4) 

where H d is the Hubble parameter at the decay and T a is the decay rate of a. 
Since H depends only on the total energy density, the decay hypersurface is 
a hypersurface of uniform total energy density, with 5Af<i = (, where ( is the 
global curvature perturbation. The equality between the sum of all energy 
densities, before the decay and after the decay, thus reads 

PA-e^ +w ^-V = p dccay = PB + e^ 1+w ^~<\ (5) 

A B 

where the subscripts — and + denote quantities defined, respectively, before 
and after the transition. In the above formula, we have used the non-linear 
energy densities of the individual fluids, which can be expressed in terms of 
their curvature perturbation (a by inverting the expression (2). 

2.1 Before the decay 

The first equality in (5) implies 

J2n A e MiA -° = 1, (3 A = 3(1 + w A ), (6) 

A 

where we have introduced the relative abundances Qa = Pa- I /Way, defined 
just before the decay. 

The above relation determines ( as a function of the (a-- Expanding it 
up to third order, one finds 



A 

A 



'Ca- + f (Ca- - C) 2 + f (Ca- - CY 



(7) 
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with the coefficients 



\a = -^-, n A = (i + WA )n A , q = J2^a. (8) 

il A 

The global perturbation ( appears on both sides of the relation (7), but one 
can use it iteratively to determine, order by order, the expression of ( in 
terms of all the ( A -, up to third order. 



2.2 After the decay 

In general, the fluid a can decay into various species, and it is convenient to 
introduce the relative branching ratio "y Aa , defined as 

7Aa = T a = ^2T Aa , (9) 

where Y Aa is the decay rate to A. Therefore, the energy density for any 
species A, just after the decay of a, is simply given by 

PA+ = PA- +lAapa- (10) 

This relation, which is fully non-linear, can be reexpressed, upon using (2), 
in the form 

e MU + -0 = (1 _ fA y A (CA-0 + f Ae Mt«-~0 ? (11) 
where the parameter 

fA = Q (12) 

represents the fraction of the fluid A that has been created by the decay. 

Expanding the above expression (11) up to third order, and using (7), 
one gets 



B 
A 

B 



c*- + ! (CB--o 2 + f (c*--cr 



-%(a + -o 2 -^(a + -o 3 , (13) 



2 V '" T " 6 

with the coefficients 



T A A = /a(i-^)a a + (1-/ a ), (14) 

T/ = /,(l-|)A. + 4 (15) 

T A C = f A (l-^\c, C^A,a. (16) 
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At second order, after substituting the first order expression for ( that 
follows from (7), Equation (13) yields 

C4+ = E t a b (b- + E (17) 



with 

A 2 



E /3eTa(5 E b - ^b){&ec - Ac) - Pa(Tab - ^b)(T ac - A c ) 



Note that the above expression for U^ c corresponds to the symmetrized 
version (with respect to the two indices B and C) of the expression given 
in [15]. This does not change the expansion (17) since U is contracted with 
a symmetric term. 

At third order, (13) yields, after using (7) up to second order, 

Ca + = E T a B Cb- + E U^Cb-Cc- + E Va CD (b-(c~( D -, (18) 

B B,C B,C,D 

with 

Vf CD = ~Y.^ eTae ^ eb - Xb ) Xf ^ 5f g- X c){5fd-Xd) 

E,F 

+ - y )^ J Pe t ae{^eb - ^b){Sec - AcX^dd - Ad) 

Ua° - - E Pe^e{$ec - ^c)(S E d - Ad) 



6 ,- 



2 

E 



—Pa(Tab — A_b) 

~~ qPa(Tab — ^b)(Tac — ^c)(Tad — Ad) • 

The above expression is the main result of this section. It provides a system- 
atic computation of the post-decay curvature perturbations for all fluids in 
a very general setting. For scenarios with several decay transitions, the per- 
turbations can be obtained by combining the various expressions of the type 
(18) for each transition. This was done explicitly, up to second order, for a 
scenario with two curvatons in [15]. In the present work, we will consider 
only single curvaton scenarios and thus use (18) only once. 

3 Mixed curvaton and inflaton scenario 

We now apply the general formalism given in the previous section to the 
mixed curvaton scenario, where fluctuations of both the curvaton and the 
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inflaton are taken into account* 2 . 

Although CDM isocurvature fluctuations and their non-Gaussianity have 
been discussed in several works [11,31-35], it is usually assumed that all of the 
CDM (or baryons for baryon isocurvature perturbations) is produced either 
before the curvaton decay or as a product of the curvaton decay. However, 
one can also envisage (as in [12,15]) intermediate situations, where the CDM 
is produced both before and during the curvaton decay. Such cases can be 
treated in the formalism discussed in the previous section and lead, as we 
will see, to interesting observational consequences. 

Here we consider a simplified scenario with three initial species: radiation 
(r), CDM (c) and a curvaton (a). We also assume that the potential for the 
curvaton is quadratic* 3 , and we thus treat a as a pressureless fluid, which 
then decays into radiation and CDM. 

3.1 Perturbations after the decay 
3.1.1 Transfer matrix 

In the present case, we find that the linear transfer matrix, whose coefficients 
are defined in (14)-(16), is given by 



where the order of the species is (r, c, a). We have used the definitions (8) 
for the coefficients A^, which leads in particular to A r = 4(1 — VL~ l ) (since 
Cl — 1 + f2 r /3 in our scenario with two nonrelativistic species). 

In the following, we will assume Q c <C 1, since the decay must occur deep 
in the radiation dominated era* 4 . As a consequence, we assume A c = and, 
therefore, x c = 0. This also implies Q = (4 — f2 CT )/3. 

In order to compare easily our results with the existing literature, it will 
be convenient to rewrite the parameter r as 



#2 The adiabatic (curvature) perturbation in such a scenario has been investigated in 



#3 The self-interacting curvaton model which include some non-quadratic potential has 
been studied in [26-30]. 

#4 Note that, although f2 c is assumed to be very small, it cannot be neglected in the 
expression for f c because -f c <j or CT can be very small, and f c can take any value between 
and 1. 




r = £r, 



(20) 



[17-21]. 
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where the quantity 



fr 



(21) 



n a 1 - (i - 7 ro .)fi C7 

can be interpreted as the efficiency of the energy transfer from the curvaton 
into radiation (£ = 1 if the curvaton decays only into radiation, i.e. 7 rcr = 1), 
and 

(22) 



is the familiar coefficient that appears in the literature on the curvaton, which 
coincides with our r only if £ = 1. 



3.1.2 Nonlinear perturbations 

It is now straightforward to obtain the post-decay curvature perturbations 
for each of the remaining fluid, up to third order, by using our general result 
(18). 

The expression for the CDM curvature perturbation reads 

C c+ = Cc+ l -US a -S c )+ l -Ul-f c ) (S a - Sc ) 2 + ±f c (l-3f c +2f c )(S a -S c )\ 

(23) 

where, on the right-hand side, we have substituted ( a — ( c = (S a — S c )/3 and 
the pre-decay subscript is implicit for all quantities. 

The above expression, which turns out to be valid also for an arbitrary 
value of x c , can in fact be derived much more rapidly by using directly (11) 
for the CDM fluid. Indeed, since (3 C = f3 a = 3, one gets 

C c+ = C + I In [(1 - / c )e 3 ^"« + f c e^~0] = Cc + ± In [l + f c (e s ^ - l)] , 

(24) 

and expanding the last expression gives immediately (23). 

The curvature perturbation for radiation is much more complicated in 
the general case. In the limit x c = 0, the full expression reduces to 



Cr+ — Cr- + 2^°"- + -^g 



21 



2r r 
3-4, + T -- 



S: 



r 

162 



9 + 18(1 - 20 1 + 4 (8£ 2 - 6£ - 3) - + 2(6£ - 1)- + 3^ 



SL (25) 



In the limit £ = 1, one recovers the expression given in [36]. 
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3.2 Primordial adiabatic and isocurvature perturba- 
tions 



In order to determine the statistical properties of the primordial adiabatic 
and isocurvature perturbations, defined at the onset of the standard cosmo- 
logical era, i.e. after the curvaton decay, one needs to relate the perturbation 
of the curvaton fluid to the fluctuations of the curvaton field, generated dur- 
ing inflation. 

For a massive curvaton without self-interaction, the relation between the 
isocurvature perturbation S a and the fluctuation 5a is simply given by the 
non-linear expression 




This result can be obtained by writing the (non-linear) energy density of the 
oscillating curvaton defined on the spatially flat hypersurfaces, characterized 
by 5N = ( r when the curvaton is still subdominant: 

p a = m 2 a 2 = m 2 (a + 5a) 2 = p^e 3 ^'^ = p a e s ° . (27) 

Expanding the expression (26) up to third order, and using the conserva- 
tion of 5a I a in a quadratic potential, we obtain 

S a ^S-\s 2 + ^S", (28) 

where the quantity 

S = 2 5 -^ (29) 

is Gaussian. 

For simplicity, we now restrict our analysis to the situation where 

Cc-=Cr- = Cinf, (30) 

by assuming that the CDM and radiation perturbations, before the curvaton 
decay, depend only on the inflaton fluctuations. This means that, effectively, 
there are only two independent degrees of freedom from the inflationary 
epoch, Ci n f and S. 

Substituting (28) into (25) and (23) then yields the primordial curvature 
perturbation 

Cr = (in i + ZlS+^Z 2 S 2 + ^Z 3 S 3 , (31) 
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with 




and the primordial isocurvature perturbation 
with 

si = f c - r, s 2 ~- 

These expressions will be useful to determine the power spectrum, the bis- 
pectrum and the trispectrum, which can be probed by observations. 



S c = Sl S+ l -s 2 S 2 + l -szS\ 



(34) 



1 / 2r 3 4r 2 

i(3/ c (l-2/ c ) + 1,-^+8^ 



r 
18 



6r 3 24r 2 
Y + l 2 ~ 



3r 



48r 
T 



15r 18 
— + 64r + — 
4 2 £ 



-36 



3.3 Power spectrum 

The power spectrum for the total curvature perturbation is given, according 
to (31), by 

= n inf + ^5 = a + A)n inf = s ~ l! ^> (37) 

where A is the ratio between the curvaton and inflaton contributions, 

A*<« (38) 

' Cinf 

and S = 1/(1 + A -1 ) represents the fraction of the power spectrum due to the 
curvaton contribution. Since the contribution from the inflaton (for standard 
slow-roll single field inflation) can be written as 

n inf = Y e Vs t" (39) 

with Vs<f, being the power spectrum of the inflaton fluctuations 5(f) and e = 
(1 /2)(dV I d(f)) 2 /V 2 being the slow-roll parameter. The parameter A is explic- 
itly given by 

8r 2 e 

A= W - (40) 
9 



Thus A depends on the curvaton model parameters, as well as on the inflation 
model. 

We can also write down the power spectrum for the isocurvature fluctu- 
ations. Using Eq. (34), we find 

V Sc = (f c - rfVs. (41) 

In our model, both curvature and isocurvature perturbations depend on the 
curvaton fluctuations. Therefore, the two types of perturbations are corre- 
lated, as quantified by the correlation coefficient: 

C = ^==e f ^ e/ = 8gn(/ c -r). (42) 

In the pure curvaton limit (S ~ 1), adiabatic and isocurvature perturbations 
are either fully correlated, if Sf > 0, or fully anti-correlated, if £/ < 0. In 
the opposite limit (S C 1), the correlation vanishes. For intermediate values 
of S, the correlation is only partial, as can also be obtained in multifield 
inflation [37]. 

Finally, it is convenient to introduce the ratio of the isocurvature power 
spectrum with respect to the adiabatic one 

This ratio can be strongly constrained by cosmological observations, but 
the precise limits depend on the assumed level of correlation between the 
isocurvature and adiabatic perturbations (since the impact of isocurvature 
perturbations on the observable power spectrum depends crucially on this 
correlation, as illustrated in [38]). Constraints on a from the latest data have 
been published only for the uncorrelated and fully correlated cases. In terms 
of the parameter a = a/(l + a), the limits (based on WMAP+BAO+SN 
data) given in [2] are* 5 

a < 0.064 (95%CL), a x < 0.0037 (95%CL) , (44) 

respectively for the uncorrelated case (S = 0) and for the fully correlated 
case (S = 1). 

According to the expression (43), the observational constraint a < 1 can 
be satisfied if \f c — r\ «r (which includes the case f c — 1 with r ~ 1) or 
if S <C 1, i.e. the curvaton contribution to the observed power spectrum is 
very small. 



* 5 Our notations differ from those of [2]: our a corresponds to their a and our fully 
correlated limit corresponds to their fully anti-correlated limit, because their definition of 
the correlation has the opposite sign (see also [39] for a more detailed discussion) . 
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3.4 Bispectrum 

We now discuss the three-point functions for our curvature and isocurvature 
perturbations, summarizing the results obtained in [15]. 

Let us start by considering, in the general case of an arbitrary number of 
observable quantities X 1 , all possible bispectra, defined by 

(J&X£Xj£> = (2tt) 3 5(£^)5 /JX (A:i, k 2 , k 3 ) • (45) 
We then assume that the X 1 can be written, up to second order, in the form 

X 1 = Nir + ^N&W + ... (46) 

(with implicit summation over the indices a and b), where the N[ and N^ b 
are arbitrary background-dependent coefficients and the 5<f) a are Gaussian 
fluctuations, generated during inflation and characterized by their two-point 
correlation functions 

(5(j) a (k)5(f) b (k')) = (2tt) 3 P ab {k) 5{k + fe') . (47) 

Substituting the decomposition (46) into the left hand side, and using (47), 
one finds [40] 

B IJK (k u k 2 ,k 3 ) = NlNiN«P a \k{)P b \k 2 ) + NlNiN«P a \k{)P^ 

+N I ab NiN«P ac (k 2 )P b \h). (48) 

In our particular case, we have only two observables X 1 = {( r , S c }, cor- 
responding to two indices which we denote / = {(, S}. Moreover, since there 
is only one Gaussian degree of freedom, S, in the non-linear expansions for 
( r and S c , respectively (31) and (34), the bispectra reduce to 

B IJK (k l , k 2 , k 3 ) = b T / L K Ps(k 2 )P § (k 3 ) + bj* 1 r Ps(h)P § (k 3 ) + bSFPgihjPgfa), 

with 

= N[ 2) N^N^, (50) 

where TVL = z 2 , = s 2 , Ni, = z±, = si, respectively. 

In order to compare these coefficients with the usual parameter* 6 /jvl 
defined in the purely adiabatic case, one must remember that f^L is propor- 
tional to the bispectrum of ( divided by the square of the power spectrum. 

#6 A11 the nonlinear coefficients that we introduce are of the local type and we drop the 
superscript "local" for simplicity. 
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The analogs of /atl can therefore be defined by dividing the coefficient b N J L K 
by the square of the ratio P^/P§ = zj/E: 

7I,JK _ 6 ,I,JK _ £^ ,I,JK ( r-,\ 
J NL — g J NL — Z A NL ■ \ 01 ) 

Taking into account the fact that the last two indices can be permuted, this 
leads to six different coefficients, explicitly given by the expressions 

fCCC _ fi^2 K,(S _ S 1^2 „ 2 ?5,CC_£^-2 ( r \ 

J NL ~ 2 " ' J NL — 3 " ' ^VL — 2 " ' v dz J 
Z l Z l Z x 

2 2 

K,SS _ s l z 2 ^ 2 ?5,C5 _ £l£2 ^ 2 /5,55 _ s l s 2 ^ 2 /rn\ 

JjVL _ 4 " ' JJVL ~~ 3 " > iVL — 4 " ' v°°/ 

z l z l z l 

It is worth noting that all the coefficients are related via the two rules 

f I,JS _ fl ,J,J< f S,IJ _ £2 f C,/J / r.A\ 

J NL — „ J NL i J NL — „ JNL ■ \°*) 

Z\ Z2 

Therefore, the hierarchy between the parameters can be deduced from the 
value of the ratios S\jz\ and s 2 /z 2 , which are given in the small r limit 
(assuming £ = 1) by the simple expressions 

- = a/^=3(^-iV ^ 3/ c (l - 2/ c ) _ 3 (r<<u = 1) . 

(55) 

Observational constraints on isocurvature non-Gaussianities are given in 
[13], for an isocurvature perturbation of the form S = S + f^S 2 , where 
S is Gaussian. The relations between the non-linear parameter and 
the parameters defined above are the following: f^ s = 2f^a> 2 , f^ L s = 
2f^ f s l > a 3 ^ 2 \C\ and = where a and C are respectively defined 

in (43) and (42). 

3.5 Trispectrum 

We now turn to the novel part of this work, which consists of the study of 
the trispectrum. 

Let us first consider the general case of several observables X 1 and let us 
introduce the trispectra T IJKL defined from the connected four-point corre- 
lation functions as 

(X^XiXgX^c = (2n) 3 5(i: i k t )T IJKL (k 1 , k 2 , k 3 , k 4 ) . (56) 
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If the observables can be written, up to third order, in the form 

X 1 = Nl8<F + 1^,5^ + + • • • (57) 

one finds that the trispectra are given by [41] 

T IJKL (k 1: fe 2 , fc 3 , fc 4 ) = Ka 2 a 3 N b J N* N^P a ' b (k 2 )P a2C (k 3 )P a3d (k 4 ) + 3 perms 

+P blC (k 3 )P aid (h)P a2b2 (k u )] + 5 perms , (58) 

with the notation fc 13 = |fei + fc 3 |, etc. 

In our particular case, where there is only one Gaussian degree of freedom 
at the non-linear level, the trispectra reduce to 

T IJKL {k u k 2 , k 3 , fc 4 ) = t^ J L KL P^k 2 )Ps{h)Ps{k A ) + 3 perms 

+ ?nl L [Ps(h)Ps(^) P s( k i 3 ) + P§(h)P s (k A )P s (ku)] + 5 perms , 

with 

41^ = NfaNfoNfcNfa, t^ L KL = NfaNfoNfcNfa, (60) 

where iV,*L = z% and = s 3 , in analogy with the notations introduced 
previously. Taking into account the symmetries under permutations of the 
indices, one finds, for two observables (I = {(,S}), 8 different parameters 
^nl L an d 9 parameters i 1 ^ 1 . 

In order to facilitate the comparison with the parameters t/vl and g^L, 
which have been defined in the purely adiabatic case by dividing the trispec- 
trum by the cube of the power spectrum, it is convenient to rescale the above 
parameters and to introduce the coefficients 

N 1 N J N K N L 

IJ,KL _ iV (2) iV (2) iV (l) iV (l) ^ 3 , v 

T NL = ~6 " ' ^ Di J 

z l 

which depend only on the second order non-linearities, and the coefficients 
~i,Jkl _ 54 r JKL _ NfoNfoN&Nfa ^ 

which depend on the third order non-linearities. 
For the purely adiabatic parameters, we obtain 

2 



r$F = |2», (63) 
& K = 4 2 S . («) 

Z l 
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which exactly coincide with the usual tnl and (54/25) #nl- The present con- 
straints on these parameters, assuming the data do not contain any isocur- 
vature contribution, are [42] 

-7.4 < 10- 5 g NL < 8.2 (95% CL), -0.6 < lO" 4 ^ < 3.3 (95% CL) . 
Similarly for the (purely) isocurvature mode, we have 

T SS,S S = fifl 3 3, (65) 

Z l 

~S,SSS s l s 3 ^3 (aa , 

9nl = — ■ ( 66 ) 
z l 

In addition to the purely adiabatic or isocurvature nonlinear coefficients, 
we also find the following cross-correlated nonlinear coefficients: 





T tt,<s _ 

'NL 


S\Z 2 „ 3 


T cs,cc _ 

'NL ~~ 


^2^2 ,-,3 

4 1—1 ' 

Z i 


(67) 




«SF c = 


S 3 ^3 


#Nl5 = 


Sl^3 ^3 
4 " ' 

Z t 


(68) 


'NL 


S 1 Z 2 ^3 
z\ " ' 


T <s,cs _ 

'NL ~~ 


S 1^2-S2 „3 

4 " 


SS,(( _ S 2 ^3 
' 'NL ~~ „4 " ' 


(69) 




s \ z ?> ^3 


sSSCC _ 
"nl ~~ 


Sl«3 ^,3 
4 " ' 




(70) 


cs,ss 

'NL 


2 

S 1^2S2 ^3 

z\ ~ 


SS,(S s l s 2 
' T NL — K 
z l 


^3 

1 — 1 5 


(71) 


~t,sss 
9nl 


S 1 Z 3 ^3 

~ ' 


p,S,(SS _ 
i/NL ~~ 


" ' 




(72) 



These nonlinear coefficients can be related by very simple rules, in analogy 
with the rules obtained for the bispectrum coefficients. The first rule is that 
the replacement of an index ( by S on the right hand side of the comma 
corresponds to a rescaling by S\j Z\. 

I.JKS _ Sl I,JKC IJ,KS _ «1 IJ,KC (7 o\ 

9NL — „ 9NL > T NL — „ NL ■ \'°) 

Zi Z\ 

The replacement of ( by S on the left hand side of the comma leads to a 
rescaling by s 2 /z2 for the t^ J £ KL and a rescaling by s 3 /z 3 for the g N J L KL - 

S,JKL _ S 3 C,JKL IS,KL_S2 K,KL /-^ 

WL — ~9NL i NL — ~~~ NL ■ 

z 3 z 2 
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These rules will be extremely useful in the subsequent discussion, enabling us 
to establish hierarchies between the 17 nonlinear coefficients defined above. 

Moreover, there are consistency relations between the trispectrum pa- 
rameters Tn£ KL and the bispectrum parameters f^^- For example, one 
finds 

1 NL — " J NL J NL i \'°) 

relating the trispectrum parameters to the bispectrum ones, which generalizes 
the purely adiabatic consistency relation tnl = (6/atl/5) 2 /5, discussed in 
[16]. 

3.6 Quantitative discussion 

In the pure curvaton scenario, when CDM is created before the curvaton 
decay (i.e. f c — 0), large isocurvature fluctuations (correlated with the adia- 
batic ones) are generated and they turn out to be too large to be compatible 
with current data. One way out is to consider the mixed inflaton and curva- 
ton scenario, where the inflaton fluctuations also contribute to the observed 
power spectrum, in addition to the curvaton ones. The isocurvature fluctua- 
tions are then "diluted," (by the S factor) and can thus be made consistent 
with observations, as studied in [11,31,34,35]. 

At the same time, since non-Gaussianity in these models originates from 
the curvaton sector, one can naively expect that non-Gaussianity will also 
be small if the curvaton contribution to the power spectrum is small. Thus 
it seems that large non-Gaussianity is difficult to realize without conflicting 
with the isocurvature constraint, although /nl > O(10) is still possible while 
satisfying isocurvature constraint [34,35]. However, in the models investi- 
gated in the present work, we also include the possibility for CDM to be 
created both from the curvaton decay and from some pre-decay epoch. As 
we will see explicitly below, this leads to very interesting consequences in the 
parameter space. 

3.6.1 Bispectrum parameters 

In Fig. 1, we have plotted the contours of and / N l in the A-r parameter 
plane, for two values of f c and assuming £ = 1. 

We recall that the parameter A, defined in (40) as the ratio between 
the curvaton and inflaton contributions to the power spectrum, is directly 
related to the parameter 5, which we have prefered to use in the analytical 
expressions for the nonlinear coefficients, by the relation 

A = j-fg. (76) 
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In the limit A 1, one gets H ~ A, whereas H ~ 1 — A -1 in the limit A ^> 1. 

The constraint (do < 0.064 at 95 % CL. ) on uncorrected isocurvature 
mode from WMAP7 is also used to identify regions still allowed by the current 
data* 7 . For small values of r, the purely adiabatic parameter is given by 



J NL 



1^2 

2r " ' 



(77) 



and the other f^ J L K can be deduced from it by appropriate factors of S\/z\ 



or S2/Z2, according to (54). 

In Fig. 2, we have plotted all the parameters f^ K as a function of f c , for 
fixed (and small) values of H and r. The figure shows how the parameters 
evolve from the region / c < r to the region f c ^> r. The figure clearly il- 
lustrates the various hierarchies between the non-linearity coefficients, which 
we discuss below. 
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Figure 1: Contours for the bispectrum coefficients and f^f S in the A-r 
plane in the cases with f c = (left) and 10~ 4 (right). The regions a > 0.064 
are shaded. 



# 7 Although this limit, strictly speaking, applies only for the region A <C 1, we have used 
the same limit for the intermediate region A ~ 1 as well as for the region A > 1. In 
the latter region, we know that the real limit is more stringent when sj = 1 since it is 
given by the constraint on a\. But constraints from the current data on a_i or as with 
intermediate values for S are not given in the literature, so we have prefered to use the 
limit on ciq everywhere as indicative constraints on the parameter space. 



16 




Figure 2: Plots of the coefficients f^ K as functions of f c . The indices I, J, K 
are specified in the figure for each curve. The other parameters are £ = 1, 
A = 10~ 3 ~ E and r = 10~ 5 . 



In the case of f c = 0, illustrated in the first plot of Fig. 1, both factors 
S\jz\ and S2/22 reduce to (—3). As a consequence, the f^ J L are slightly 
enhanced with respect to by a factor (— 3) /s where Is is the number 
of indices S. In particular the pure isocurvature parameter is (—27) times 
Jnl as confirmed by the figures. The constraint on the isocurvature power 
spectrum imposes 5 1 and the non-Gaussianity can therefore be significant 
only if r is sufficiently small to compensate the E 2 suppression. 

When f c does not vanish, as illustrated in the second plot of Fig. 1, a new 
region appears in the parameter space, where the isocurvature constraint is 
satisfied even if E ~ 1. This requires a fine-tuning between f c and r, more 
precisely, that of 

f c -r^e f y^r (S~l), (78) 

according to (43). With respect to the purely adiabatic non-Gaussianity, 
the other types of non-Gaussianities are suppressed in general (except f N 'i , 
which could be of the same order of magnitude in the particular case £ 1 
and f = r/£~l). 
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In the case r <C f c <C 1, another interesting feature appears in the region 

S SS 

S< 1. Indeed, the isocurvature bispectrum amplitude / NL becomes larger 
than /nlS as pointed out in [15]. The reason is that the factors S\j Z\ and 
s 2 /z2 become 3/ c /r in this region so that the J'nl with adiabatic indices are 
suppressed with respect to the purely isocurvature non-Gaussianity. In this 

C CC S SS 

region, and / NL are respectively given by 

K,CC ^ _ J^ a 2 is,ss _ ^/?~2 _ ^ (r f 1} (791 

/nl ~2r" "54/ c 4 ' /nl - 2r 4" - 2/c ' {r^j c ^l),{(V) 

from which we can see that f^f S is enhanced by the factor (3/ c /r) 3 compared 
to /nl and can be significant if / c is well below a 2 . In such a scenario, future 
observations would thus detect non-Gaussianity from the isocurvature per- 
turbations rather than from the adiabatic ones. For the correlated nonlinear 
coefficients f^ K , we obtain 

fCXS „ fS,a „ ^fc^2 K,SS „ fSXS „ ^1H^2 (m) 

J NL J NL — 2 r 2" ' ■'NL ~ J — 2 r 3 " ' 

Their value depends only on the total number of adiabatic indices. The 
hierarchy between the non-linearity parameters depends only on the ratio 
fc/r. 

In the opposite limit where f c is much smaller than r (/ c < r <C 1), 
all the non-linearity coefficients become independent of f c . Once again, the 
amplitude of the purely isocurvature coefficient is enhanced (now by a factor 
27) with respect to the purely adiabatic one, but the sign is changed: 

/nl - 2r ^ - 54r , /nl - 2r ^ • 
The correlated nonlinear coefficients have intermediate values, 

K,CS „ fS,C( _ _$_^2 K,SS _ rS,(S _ 27 „ 2 ( n 2) 

J NL — J NL — ' ''NL — ^NL — ^ " ■ l oz J 

All these values differ simply by powers of (—3), since Si/zi ~ 52/^2 — — 3 in 
this limit. An interesting consequence is that the nonlinear coefficients with 
an odd number of isocurvature indices are negative. 

3.6.2 Trispectrum parameters 

The contours of the purely adiabatic and isocurvature nonlinear coefficients 

CC CC SS SS C CCC 

and r NL ' have been plotted in Fig. 3, while the coefficients <7nl an d 
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Figure 3: Contours for the trispectrum coefficients and t^'1 ,ss in the 

A-r plane for the cases f c = (left) and IO -4 (right). The regions a > 0.064 
are shaded. 



9nl SS are depicted in Fig. 4. In these figures, the region excluded by current 
observations of the isocurvature mode is again shaded. 

For small values of r, the purely adiabatic trispectrum coefficients are 
(for £ = 1) 

4^ * ^H 3 , ~ -^~ 3 , (r « 1, £ = 1). (83) 

They can be also related to the purely adiabatic nonlinear coefficient J'nl = 
fNX by the following expressions [16]: 

T NL -JNlI^ 9nL — 7^JNL^- (?V 

This implies that, when S C 1, ^vx * s enhanced with respect to jfjyx 
whereas, by contrast, 9nl ls suppressed with respect to /jvl- 

For / c = 0, all Tj/l are slightly enhanced with respect to the adiabatic 
r NL by factors (— 3) /s . In particular, t^ S £ SS = 81t nl . The same conclusion 
applies to the g 1 ^^ 1 since S3/Z3 = —3 for f c = 0. 

In the fine-tuned region f c ~ r, the purely isocurvature nonlinear coeffi- 
cients and most of the correlated coefficients are suppressed with respect to 
the purely adiabatic ones (except in the case £ -C 1, as discussed earlier). 
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Figure 4: Contours for the trispectrum coefficients g^L** and g^f SS in the 
A-r plane for the cases with f c = (left) and 10~ 4 (right). The regions 
a > 0.064 are shaded. 



However, since s^/zs ~ 3/2, one finds 

//vr - I'lXr- (/e~r, 5-1, $ = 1). (85) 

In the region r < / c < 1 and 5<C 1, where the bispectrum is dominated 
by its purely isocurvature component, we find the same conclusion for the 
trispectrum coefficients since S3/Z3 ~ f 2 /r 2 . The largest coefficients are thus 
(for £ = 1) 

ssss 729 f 4 a 3 s sss 2187/* , 3a 3 „ 

nl — A « - — A tii 9nl - „fi - — « (r «^ j c «^ i,<^ - rj, 



4r 6 ~ 4/ c 2 ' yJVL ~ 2r 6 ~ ~ 2/ c 

SSS 

and one finds the following relations with f N ' L : 



(86) 



t nl - yf nl ) a > ^Wl - -3a/jvt , (87) 

which are very similar to their adiabatic counterparts (84). 

In Fig. 5, we have plotted the evolution of the hierarchy between all the 
parameters t^ KL as the parameter f c varies, r being kept fixed (with a small 
value). In the limit r < / c < 1, one sees clearly that the purely isocurvature 
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Figure 5: Plots of the coefficients t^ kl as functions of f c . The indices 
I, J, K, L are specified in the figure for each curve. The other parameters are 
f = 1, A = 10" 3 and r = 1(T 5 . 



coefficient dominates, while the other coefficients vary according to 

4t KL * (j^ k 4 S £ SS , (r « fc « 1), (88) 

where 1^ is the number of adiabatic indices among the four indices of the 
coefficient. In this opposite limit, / c < r < 1, we have 

ij,kl _ , / 1 \ ss,ss 



4t .-{-3) t nl, (/c«r«l). (89) 

In Fig. 6, we have plotted the evolution of all the g^ KL as f c varies. Once 
again, the two opposite limits have very distinctive hierarchies between the 
nonlinear coefficients. In the limit r C / c C 1, 9nl SS dominates. Moreover, 



§&r--^s^g^r, (90) 



and all the other coefficients can be deduced by using the relations 

~S,JKL ( r \ C ~S,SSS -CJKL ( r \ C ~(,SSS / f 

9nl - ( 3J ) 9nl , 9nl - ( gj ) > ( r < /c < 1), 

(91) 
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where /| is the number of adiabatic indices among the three indices after the 
comma. In the opposite limit, / c C r C 1, the hierarchy is simply 

~g\ii KL * (--) k ?nl 8S > (/c « r « ^ ( 92 ) 

where is the number of adiabatic indices among the four indices. 



4 Conclusion 

In this paper, we have investigated non-Gaussianity in models with isocur- 
vature fluctuations, paying particular attention to the trispectrum. After 
presenting a general formalism to calculate bi- and tri-spectra of adiabatic, 
isocurvature and correlated types, allowing various decay scenarios, we have 
applied the formalism to the mixed curvaton and inflaton model. 

We have studied how the amplitude of the non-linearity parameters, 
which consist of six f^ K , nine t^ kl and eight L , depend on the param- 
eters r, the curvaton contribution to the power spectrum S (or, equivalently, 
A) and f c , which represents the fraction of dark matter produced from the 
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curvaton decay. We have found that large non-Gaussianity (in the bispec- 
trum and trispectrum) from both of adiabatic and isocurvature modes is 
possible in some cases without conflicting the isocurvature constraint from 
the measured CMB power spectrum. We have also compared the relative 
size of the non-linearity coefficients of various types as depicted in Figs. 2, 5 
and 6, and found that different regions in the parameter space correspond to 
very distinctive hierarchies between the non-Gaussianity coefficients. 

Observations of non-Gaussianity, and in particular of the trispectrum 
have become an important goal in cosmology. From another perspective, 
isocurvature fluctuations are associated with the generation of dark matter 
and baryon asymmetry in the Universe. Thus, non-Gaussianity from isocur- 
vature fluctuations, if detected in the future, would give us a lot of insight 
into the nature of dark matter, the mechanism of baryogenesis, and therefore 
into high energy physics. 
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